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Light quark masses for dynamical, non-perturbatively 0(a) improved 
Wilson fermions 

D. Plciter a ' b , QCDSF and UKQCD Collaboration 

a Deutsches Elektronen-Synchrotron DESY & NIC, D-15735 Zeuthen, Germany 
b Institut fur Theoretische Physik, Freie Universitat Berlin, D-14195 Berlin, Germany 

We present results for the light quark masses in lattice QCD with two degenerate flavours of dynamical fermions. 
We used configurations generated by the UKQCD and QCDSF collaborations at six different combinations of /3 
and K sea . 



1. INTRODUCTION 

The calculation of the light quark masses is one 
of the major aims of lattice gauge theories. In 
recent years these fundamental parameters have 
been determined with a rather high accuracy us- 
ing quenched QCD. At least for the strange quark 
mass, systematic errors are believed to be un- 
der control, except for quenching effects. How- 
ever, for current simulations of QCD with dy- 
namical fermions this is not yet the case. The 
main sources of uncertainty are the lack of non- 
perturbative results for the renormalization con- 
stants and the missing possibility to do a reliable 
continuum extrapolation. 

For the calculations presented here we used dy- 
namical configurations with Nf = 2 flavours of 
degenerate quarks for six different combinations 
of j3 and K sca . The coefficient for the improve- 
ment term of the Wilson fermion action, csw(fl), 
was taken from The lattice spacing a varies 
between 0.104 and 0.089 fm. For the lightest sea 
quarks the ratio mps/wv was 0.6. (For further 
details, see |0j.) 

We define the renormalized quark mass in the 
MS-scheme at scale fj, in the following way: 
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where 



MS 
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(l + b A am q )Zf s 
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and 



1 (G> 4 A 4 ft)P(0)) 

2 (P(t)P(O)) 



(3) 



are the bare quark masses from the conserved vec- 
tor current (CVC ) and the partially conserved 
conserved axial vector current (PC AC ), respec- 
tively. In order to eliminate 0(a) effects intro- 
duced by the axial vector current A„, one has to 
replace A 4 in the last equation by 



A 4 — > A4 + caclO^P, 



(4) 



where P is the pseudoscalar density. We there- 
fore have to compute the following ratios on the 
lattice: 



(5) 



(d 4 A 4 (t)P(0)) 



ca 



(adlA 4 (t)P(0)) 



2(P(t)P(0)) " 2{P(t)P(0)) ■ 

Since both ratios can be measured rather pre- 
cisely, we are, for the rest of this talk, left with 
basically two problems: 

• Determination of the renormalization con- 
stants Zq, the coefficients bo, which 
parametrize the quark mass dependence of 
the renormalization, and the improvement 
coefficient ca- 

• Extrapolation of the lattice results to the 
physical quark masses. 

In this talk we use ro = 0.5 fm and take the 
lattice results for r^/a from M to set the scale. 
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2. IMPROVEMENT AND RENORMAL- 
IZATION 

Until now, none of the renormalization con- 
stants and improvement coefficient, except for 
Za 0], have been calculated non-perturbatively 
for N{ = 2 improved fermions. We there- 
fore will use perturbative results and compare 
these with available non-perturbative numbers, 
e.g. quenched results, in order to have some esti- 
mate of the possible errors. 

Ordinary perturbation theory on the lattice is 
known to suffer from problems because of large 
tadpole diagrams and the expansion in a non- 
physical (bare) coupling constant. Usually it is 
better to apply so-called tadpole improved per- 
turbation theory, which refers to a strategy to 
sum tadpole diagrams and to use a physical cou- 
pling constant in the perturbative expansion. 

The choice of the coupling constant to be used 
is not unique. We will use the boosted coupling 
constant g* 2 = g 2 ju\ for bare quantities and the 
coupling constant in the MS-scheme (g MS ) 2 = 
Air a^ s (f/a) for the computation of renormal- 
ized quantities. a^ s (l/a) can be calculated from 
the 4-loop expansion of the p- function Q , using 
A m = 250 MeV. 



convenient to define the ratio 



;(!) 



Ri(a) 



+ ca m y q 



(2) 



4 NP) 4 2) 
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In Fig. Jfi we show Ri(a) using our quenched 
results Mf. In the continuum limit we expect 
i?i(0) = I. Although R\(a) approaches this 
limit very quickly, the results for the bare quark 
mass show discretization effects of up to 15% for 
a ~ 0.1 fm. 
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Figure 2. Chiral extrapolation of m q as func- 
tion of 1/K val for constant (/3, K SGa ). The non- 
perturbative and the tadpole-improved results for 
k™ 1 are marked with x and o, respectively. 



R,(a) 



(a/r ) 



Figure 1. Ratio R\(a) defined in Eq. (^|) as a 
function of the lattice spacing using ca = (O) 
and the tadpole- improved result (♦). 



The improvement coefficient ca is known from 
1-loop perturbation theory ||. In the quenched 
approximation ca has also been computed non- 
perturbatively To parametrize the difference 
between the non-perturbative and the tadpole- 
improved results at the strange quark mass, it is 



As a next example we compare the tadpole- 
improved result for k™ 1 , which is given to 1-loop 
tadpole-improved perturbation theory by Q 

■f = \ [1 + g* 2 (0.025 - 0.029 c sw u 3 

8 (7) 

-0.012 (c S wUo) 2 )K\ 



with the non-perturbative results. For each 
K sea ) we calculated m q for three different val- 
ues of K val . From a linear extrapolation in 1/ re val 
to m™ 1 = 0we determine k™ 1 . In Fig. |^ we com- 
pare our non-perturbative results with the results 
from Eq. (^). In this case we find the relative dif- 
ference between both results to be less than 0.1%. 

The renormalization constants are known to 
1-loop tadpole-improved perturbation theory ||. 
Z^ s (fi) is scale dependent and we will use the 
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Figure 3. Non-perturbative (♦) and tadpole- 
improved results (O) for Zp/ZsZa as a function 
of the pseudoscalar mass at K sei 



„ val 



renormalization group equations, to scale our re- 
sults to a convenient scale, i.e. 



•MS 



(2 GeV) = 



AZ MS (l/a) ™ 



AZ MS (2GeV) 



fnf (1/a), (8) 



where AZQjl) is the factor which translates 
a scale and scheme dependent value of the quark 
mass to its renormalization group invariant value. 
To first order we expect that the following re- 



lation between m q and m q holds: 



Zv 



ZsZa 



0{ml) 



(9) 



From the slope of the fits shown in Fig. g we get 
a non-perturbative value for the ratio Zp / Z§Za- 
In Fig. H we compare these results with the cor- 
responding tadpole-improved numbers. In this 
case we find the difference to be of the order of 
10%. However, one might hope that the results 
for Z m (ji) — Za/Zp(/j,) from tadpole- improved 
perturbation theory to be more reliable, since 
the contributions from the tadpole diagrams can- 
cel. For the quenched approximation this can be 
checked, since Z m (fi) has been calculated non- 
perturbatively in the Schrodinger functional (SF) 
scheme [ fTo| . In order to compare perturbative 
and non-perturbative results we looked at the ra- 
tio 

AZ M (l/fl)Z,f(l/ fl ) 



R 2 (a) 



AZ**(l/L)Zg(l/L) 



(10) 



where we expect i?2(a) = 1 + 0(a). The results 
for i?2(«) are shown in Fig. [|. 

In perturbation theory we have 6a ~ bp- This 
has recently been confirmed non-perturbatively 



for quenched QCD |[Ll| . It is therefore expected 
that these coefficients will have only a small effect 
on the results. 
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Figure 4. Ratio i?2(a) defined in Eq. (|Io| ) for 
quenched QCD at /3 = 6.0,6.2, and 6.4 using 
results from ordinary (solid line) and tadpole- 
improved perturbation theory (♦). 



3. QUARK MASSES 

We are now able to combine our results and 
to calculate m^ s (2GeV). In Fig. | we show the 
squared pseudoscalar mass as a function of the 
renormalized quark mass at K val = re soa . 

Even for large quark masses the deviation from 
linearity is very small. Making a linear ansatz we 
can compute the light and strange quark mass by 
solving the equations 



bi 2r$m\ 
h {r fh\ + r fh s ) 



(ro< hys ) 2 , 
(r m£- hys ) 2 . 



(11) 
(12) 



Using the experimental 
137.3 MeV and 



values mP hys 



-.phys 



495.7 MeV yields 
3.8(3) MeV and 



the results m 1 Mb (2GeV) 
m;p(2 GeV) = 96(4) MeV. 

For a more sophisticated ansatz we make use of 
partially quenched chiral perturbation theory. To 
1-loop order the pseudoscalar mass can be written 
in the following form |L2] : 



m 



PS 



6!( K soa ) 2rn; al + b 2 ( K sca ) (2m™ 1 ) ' 



(13) 



The coefficients b\ and bi can be calculated from 
a quadratic fit to our results at constant (/3, K sca ). 
The results, which show a rather mild dependency 
on K sca , are plotted in Fig. [|. To calculate the 
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2r m q (2GeV) 
Figure 5. Renormalized quark mass as function 
of the squared pseudoscalar mass in units of ro . 

light and strange quark masses we use the ansatz 

(r mP^ s ) 2 := &x« a ) r 2m + 6 2 (C) fa 2mi) 2 
(r r< hys ) 2 :=6 x ( K f a ) r (m 1 + m s ) + 
& 2 « a ) [r (mi + m s )] 2 . 

We finally get the following results: 

mp(2GeV) = 3.5(2) MeV, 



(14) 



15) 



-MS 



(2GeV) = 90(5) MeV. 



m 



In a previous quenched calculation JTj we found 
(2GeV) = 105(4) MeV. Since systematic er- 



MS 



rors, which are currently not under control, might 
be of the order 10%, we can not conclude that 
the (small) difference between these results is a 
quenching effect. 
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Figurc 6. Coefficients b± and 62 as a function 
of the squared pseudoscalar mass at K sca = n . 
The lines show a linear fit. 



4. CONCLUSIONS 

We did a preliminary analysis of the light 
quark masses using two flavours of degenerate, 
non-perturbatively improved fermions. We used 
tadpole-improved perturbation theory, but tried 
to give some idea of the systematic errors. For the 
extrapolation to the physical region we used par- 
tially quenched chiral perturbation theory. The 
difference between our final results and earlier 
published quenched results was found to be small. 
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